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Abstract 

Investigations on hippocampal and amygdalar volume have revealed inconsistent results in patients with posttraumatic 
stress disorder (PTSD). Little is known about the structural covariance alterations between the hippocampus and amygdala 
in PTSD. In this study, we evaluated the alteration in the hippocampal and amygdalar volume and their structural covariance 
in the coal mine gas explosion related PTSD. High resolution T1 -weighted magnetic resonance imaging (MRI) was 
performed on coal mine gas explosion related PTSD male patients (n = 14) and non-traumatized coalminers without PTSD 
(n = 25). The voxel-based morphometry (VBM) method was used to test the inter-group differences in hippocampal and 
amygdalar volume as well as the inter-group differences in structural covariance between the ipsilateral hippocampus and 
amygdala. PTSD patients exhibited decreased gray matter volume (GMV) in the bilateral hippocampi compared to controls 
(p<0.05, FDR corrected). GMV covariances between the ipsilateral hippocampus and amygdala were significantly reduced in 
PTSD patients compared with controls (p<0.05, FDR corrected). The coalminers with gas explosion related PTSD had 
decreased hippocampal volume and structural covariance with the ipsilateral amygdala, suggesting that the structural 
impairment of the hippocampus may implicate in the pathophysiology of PTSD. 
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Introduction 

Posttraumatic stress disorder (PTSD) is an anxiety disorder 
observed in people who have been exposed to severe emotional or 
physical life-threatening traumatic events (such as war, sexual 
abuse, and natural disasters), and reacts with intense fear, 
helplessness or horror. PTSD is characterized by multiple chronic 
and disabling symptoms, including re-experiencing the event, 
avoidance of clues related to the trauma, and hyperarousal. It has 
been proposed that the structural and functional changes of the 
emotional neurocircuitry, mainly consisted of the amygdala, 
medial prefrontal cortex and hippocampus, may contribute to 
the development of PTSD [1]. 

The relationship between the amygdala and the medial 
prefrontal cortex was investigated deeply in PTSD by many 
consistent studies, demonstrating that one candidate neural 
mechanism of PTSD may be the failure of medial prefrontal 
cortex in suppressing the hyperresponsivity of the amygdala to 
threat-related stimuli, which mediates a core symptom of 
hyperarousal of PTSD [1,2]. The PTSD is also reported to be 
closely related to the dysfunction of emotional memory processing. 
The amygdala plays a critical role in memory of stressful events 
[3], and the hippocampus is involved in affective processing [4,5]. 
The structural alteration of hippocampus has received abroad 
attention in PTSD. However, the results on gray matter volume 



(GMV) of the hippocampus in PTSD are controversial. Although 
many studies reported a reduced hippocampal volume in PTSD 
patients compared with Non-PTSD subjects [6-14], several 
studies failed to find such GMV differences [15-18]. Several 
meta-analyses confirmed reduced hippocampal volume in PTSD 
patients; however, the laterality effects are inconsistent. For 
example, some studies reported GMV reduction in the bilateral 
hippocampi [6,19,20], and others revealed GMV reduction only 
in the right hippocampus [21], or in the left hippocampus [22], in 
PTSD patients. The heterogeneous properties of the PTSD 
subjects, including gender, age, trauma type, duration of trauma 
exposure, severity of trauma, and elapsed time since trauma 
[6,19,23], may account for the aforementioned inconsistent results 
regarding hippocampal volume in PTSD. The reports on the 
structural changes in amygdala in PTSD were fewer and also 
inconsistent. Increased [24], and decreased [25,26], amygdala 
volume was both reported in PTSD, Two meta-analyses both 
showed on significant alteration of amygdala volumes in PTSD 
[27,28]. 

Relative to the structural alterations of the hippocampus and 
amygdala, the structural covariance between the two regions has 
not been reported in PTSD. The interaction between hippocam- 
pus and amygdala is important for processing the emotional 
memory under normal condition [29-31]. The amygdala and the 
hippocampus have real anatomical connectivity [32,33]. A pilot 
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study revealed that the GMV of the hippocampus covaries 
strongly with the amygdala in healthy subjects [34]. Structural 
covariance quantifies the extent to which structure covaries among 
different brain regions across individuals. Although the biological 
meaning of this structural covariance remains controversial, it 
appears to reflect mutually trophic influences or common 
experience-related plasticity [35,36], and may reflect synchronized 
development of covarying brain regions. Covariance methods can 
be theoretically applied to any structural or anatomical phenotype 
(such as surface area, GMV, and diffusion indices, etc.) and could 
provide additional information regarding the functional or 
structural networks. Networks of structural covariance partially 
recapitulate the functional networks of healthy subjects and allow 
the exploration of abnormal structural covariance networks related 
to various neurological and psychiatric diseases [36]. Therefore, 
evaluating the alteration of structural covariance between the 
amygdala and the hippocampus may broaden our insight into the 
pathogenesis of PTSD. 

The structural changes of the hippocampus and amygdala in 
PTSD patients resulted from combat, natural disaster, child abuse, 
and civilian assaults were extensively investigated. However, only 
one research group had reported the hippocampal structural 
alteration in coal mining disaster-related PTSD. They found that 
coal mine floods-related PTSD subjects had significantly reduced 
fractional anisotropy value in bilateral hippocampal body [37], 
and decreased volume and density in the left anterior hippocam- 
pus compared with survivors without PTSD [38]. The victims of 
the coal mine disaster had high homogeneity in demographic 
background, trauma intensity, and duration of trauma exposure, 
which offers an advantage in evaluating the PTSD-related 
structural brain damage. 

In this study, we recruited a group of male coalminers with coal 
mine gas explosion-induced PTSD and a group of matched non- 
traumatized coalminers to investigate the alterations in GMV of 
hippocampal and amygdala, and further investigate alterations in 
GMV covariance between the hippocampus and the amygdala in 
PTSD. We hypothesized that the hippocampal volume was 
reduced, and GMV covariance between hippocampus and the 
ipsilateral amygdala was disrupted in PTSD patients who have 
suffered from a coal mine gas explosion. 

Methods 

Subjects 

A coal mine gas explosion disaster occurred in February 2006 in 
Shanxi province of China. Thirty male survivors from the disaster 
were diagnosed as PTSD according to the Diagnostic and 
Statistical Manual of Mental Disorders (DSM-IV) criteria at 6 
months after trauma. Seven years later, 24 of the 30 victims with 
PTSD were recruited again for this MRI study. They were re- 
evaluated by a trained psychiatrist to ensure the persistence of 
PTSD. Clinician-Administered PTSD Scale (CAPS) was used to 
assess symptom severity of each patient. Twenty five male 
coalminers without exposure to the trauma, the colleagues of the 
victims, were enrolled as controls. 

Conventional MR images were used for excluding subjects with 
visible brain lesions. All controls exhibited a normal MRI 
appearance. Eight PTSD patients were excluded for visible brain 
lesions (encephalomalacia, tumors, or cysts). Two patients were 
further excluded because they met diagnostic criteria only for the 
lifetime but not for the current PTSD. Only the remaining 14 
patients with current and lifetime PTSD were used for further 
analysis. For each subject, anxiety and depression symptoms were 
assessed using the Hamilton Rating Scale for Anxiety (HAMA) 



and Depression (HAMD-24 items version), respectively. The 
experiment was approved by the Ethical Committee of Tianjin 
Medical University General Hospital and written informed 
consent was obtained from each subject before the study. 

MRI acquisition 

MR images were obtained on a 3.0 T scanner (Magnetom 
Verio, Siemens, Erlangen, Germany). A Tl -weighted volumetric 
magnetization-prepared rapid gradient-echo sequence was used to 
acquire a series of 188 contiguous sagittal high resolution 
anatomical images with the following parameters: repetition 
time = 2000 ms, echo time = 2.26 ms, inversion time = 900 ms, 
flip angle = 9, matrix = 256 x224, field of 

view = 256 mmx224 mm, slice thickness = 1 mm, which resulted 
in an isotropic voxel of 1 mm x 1 mm x 1 mm. 

MR data preprocessing 

Structural images were preprocessed and analyzed using the 
Statistical Parametric Mapping software (SPM8; http://www.fil. 
ion.ucl.ac.uk/spm/software/spm8). The structural MR images 
were segmented into gray matter (GM), white matter and 
cerebrospinal fluid. GM images were subsequendy spatially 
normalized into a Montreal Neurological Institute (MNI) space 
(http://www.mni.mcgill.ca/) with a resolution of 
1.5 mm x 1.5 mm x 1.5 mm using the diffeomorphic anatomical 
registration through exponentiated lie algebra (D ARTEL) method 
(Ashbumer, 2007). GMV of each voxel was obtained by 
multiplying the GM concentration map by the non-linear 
determinants derived during spatial normalization. Here, the 
regional GMV represents normalized GMV after removing the 
effect of variance in individual brain sizes. Finally, the GMV 
images were smoothed with a full width at half maximum 
(FWHM) kernel of 8 mm. After spatial preprocessing, the 
normalized, modulated, and smoothed GMV maps were used 
for analysis of hippocampal volume and structural covariance. 

Voxel-wise GMV comparisons in hippocampus 

The bilateral hippocampi and amygdalae were extracted using 
the Anatomical Automatic Labeling (AAL) atlas (www.cyceron.fr/ 
freeware/); and they were combined and used as the mask during 
inter-group comparisons (Figure 1). Inter-group differences in 
GMV were voxel-wise compared within the mask using a two- 
sample Hest. A false discovery rate (FDR) method was used to 
correct for multiple comparisons, and a corrected threshold of p< 
0.05 as well as a minimum of 100 contiguous voxels was 
considered statistically significant. 

Clusters with altered hippocampal GMV in PTSD patients 
were defined as regions of interest (ROIs). The GMVs in the 
hippocampal ROIs were extracted and the inter-group differences 
in these ROIs were further compared using a two-sample Hest, 
and the Cohen's d was used to represent the effect size of the 
comparison. Considered that the anxiety and depression might 
confound the result, the one-way Analyses of Covariance 
(ANCOVA) controlling for HAMD and HAMA score (both 
separately and concurrently) were further performed to compare 
the inter-group differences of GMV in the hippocampal ROIs. A 
partial correlation analysis controlling for age effect was used to 
test correlations between GMV of the hippocampal ROIs and the 
PTSD symptom severity (the total and subtype scores of CAPS). 
Significance level was set at/><0.05. 
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Figure 1. Combined hippocampal and amygdala mask on 
sagittal planes. Inter-group comparison of gray matter volume is 
performed within this mask. Left, left hemisphere; right, right 
hemisphere. The below numbers are the Montreal Neurological 
Institute coordinates on x axis. 
doi:1 0.1 371/journal.pone.01 02042.g001 

Structural covariance analysis 

Voxel-wise GMV covariance quantifies the extent to which 
GMVs covary between different brain regions across participants. 
Here, the seed regions were defined as hippocampal ROIs that 
exhibited significant inter-group differences in GMV. The 
bilateral amygdalae were extracted with AAL atlas and were used 
as the mask. The GMV covariance between each seed region and 
voxels within ipsilateral amygdalar mask was obtained by 
calculating correlation coefficients between the GMV of the seed 
region and those of the ipsilateral amygdala voxels across all 
subjects of each group. Significance level was set at ^<0.05, FDR 
correction for multiple comparisons, and with a minimum of 100 
contiguous voxels. 

A one-way ANCOVA was used to detect the inter-group 
differences in GMV covariance of the hippocampal ROIs with the 
ipsilateral amygdala (/j<0.05, FDR corrected, cluster size > 
100 voxels). ROI-based analysis was also performed to validate the 
results of the voxel-wise analysis. Here, the amygdalar ROIs were 
defined as the amygdalar clusters exhibited significant inter-group 
difference in GMV covariance with the hippocampal ROIs. GMV 
correlations between the ipsilateral hippocampal and amygdalar 
ROIs were examined using the Pearson correlation analysis (p< 
0.05) in two groups, respectively. 

Statistical analysis for demographic and clinical variables 

Kolmogorov-smirnov tests were performed to analyze the 
distribution of demographic and clinical variables. Differences 
between groups were analyzed with the Student's t test in case of 
normally distributed variable (age) or with Mann-Whitney U test 
in case of non-normally distributed variables (HAMA and HAMD 
score) in SPSS 11.0 (SPSS, Inc, Chiago.IL). Significance level was 
set at P<0.05. 

Results 

Demographic and clinical data 

All subjects recruited in this study were male coalminers who 
had the same socioeconomic state and working environment. 
There were no significant differences in age and HAMA score 
between the two groups (/)>0.05), but the HAMD score in PTSD 
patients was higher than that in controls (p<0.001). All non-PTSD 
subjects and 1 1 PTSD patients had the anxiety state (HAMA score 
>14). One non-PTSD subjects and 13 PTSD patients had the 
depression state (HAMD score —8). All the PTSD patients had no 



substance abuse before and after traumatic event. The demo- 
graphic and clinical data are shown in Table 1. 

Reduced hippocampal volume 

Coalminers with PTSD exhibited decreased GMV in the 
bilateral hippocampi compared to non-traumatized coalminers 
(Figure 2A, Table 2). ROI-based analysis showed that the 
volumes within the left and right hippocampal ROI decreased 
7.96% (controls: 0.653±0.039, PTSD: 0.601 ±0.039; /5<0.001, 
Cohen's d = 0.95) and 7.26% (controls: 0.496±0.033, PTSD: 
0.460±0.031; /> = 0.002, Cohen's d = 0.80) respectively in PTSD 
patients compared to controls (Figure 2B). The ANCOVA shows 
intergroup differences of GMV in the left hippocampal ROI was 
still significant after controlling the HAMD and HAMA score 
(both individually and concurrently). However, the intergroup 
difference of GMV in the right hippocampal ROI disappeared 
after controlling the HAMD scores or HAMD and HAMA scores 
concurrendy (Table 3). There were no significant inter-group 
differences in GMVs of the bilateral amygdalae. In PTSD patients, 
there was no any statistical correlation between the hippocampal 
volume and the HAMA, HAMD, and CAPS score. 

Structural covariance between the hippocampus and the 
ipsilateral amygdala 

Voxel-wise structural covariance analysis revealed that the 
GMV of the left hippocampal ROI significantly covaried with the 
left amygdala (MNI coordinate: -20, -9, -17; voxels: 478; t= 9.08, 
j6<0.001, FDR corrected), and that the GMV of the right 
hippocampal ROI significandy covaried with the right amygdala 
(MNI coordinate: 32, -4, -27; voxels: 452; t-AAQ, /> = 0.012, 
FDR corrected) in controls (Figure 3). However, there was no 
statistically significant GMV covariance between hippocampal 
ROIs and the ipsilateral amygdala in PTSD patients. 

Intergroup comparison of hippocampal structural covariance 
revealed a significandy reduced GMV covariance between the left 
hippocampal ROI and the left amygdala cluster (MNI coordinate: 
-21, -3, -21; voxels: 450; t= 4.87, p<0. 001, FDR corrected), and 
between the right hippocampal ROI and the right amygdala 
cluster (MNI coordinate: 33, -3, -29; voxels: 374; /=2.74, 
p — 0.042, FDR corrected) in the PTSD patients compared with 
controls (Figure 4). ROI-based correlation analysis showed strong 
positive correlations between the GMVs in the left hippocampus 
and the left amygdala (r= 0.776, /><0.001) as well as between the 
GMVs in the right hippocampus and the right amygdala 
[r= 0.568, p = 0.003) in healthy controls, while these GMV 
covariance disappeared in PTSD patients (left side: r— -0.205, 
p = 0.483; right side: r = -0.219, p = 0.453) (Figure 5). 

Discussion 

In the present study, we investigated the alterations of the 
hippocampal and amygdalar volume and structural covariance 
between the hippocampus and the amygdala in coal mine gas 
explosion-related PTSD patients. We found that coalminers with 
PTSD had a significandy decreased GMV of hippocampus and 
decreased GMV covariance between the hippocampus and 
ipsilateral amygdala compared to non-traumatized coalminers, 
which may be associated with the dysfunction of emotional 
memory processing in PTSD patients. 

Reduced hippocampal volume in PTSD patients 

Reduced volumes in bilateral hippocampi were revealed in 
coalminers with PTSD, especially in the left hippocampus, which 
is consistent with previous meta-analyses that reported a close 
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Table 1. Demographic information of the PTSD patients and controls. 







PTSD (n = 14) 


Controls (n = 25) 


statistics 


P 


Age (years) 


33.1 ±5.4 


36.0±5.7 


f= -1.497 


0.143 


HAMA 


15.5(5-23) 


15(14-20) 


z= -0.400 


0.689 


HAMD 


15(5-19) 


2(0-10) 


z= -5.059 


<0.001 


CAPS lifetime 


87.6+14.1 








CAPS current 


50.4±18.8 









HAMA and HAMD scores were non-normally distributed and showed as median (range). CAPS, Clinical-Administered PTSD scale; HAMA, Hamilton rating scale for 
anxiety; HAMD, Hamilton rating scale for Depression; PTSD, posttraumatic stress disorder. 
doi:1 0.1 371 /journal.pone.01 02042.T.001 



association between the hippocampal volume and PTSD, 
although causal relationships between them remain unclear 
[6,19,20]. A study found that combat-exposed severe PTSD 
patients and their twins both had significantly smaller hippocam- 
pal volumes than non-PTSD twin pairs, suggesting that small 
hippocampal volume may be a risk factor for developing PTSD 
following trauma [39]. On the contrary, other studies demon- 
strated that the smaller hippocampal volume was a consequence of 
PTSD. Studies of animal model confirmed that traumatic 
experience could lead to reduced hippocampal volume [40,41]. 
The reduction in hippocampal volume can be explained by the 
increased release of glucocorticoid under trauma exposure, which 
can affect hippocampal structure by causing atrophy of dendritic 
processes, neuronal loss, and decreased levels of brain-derived 
neurotrophic factor [42]. A longitudinal study on PTSD also 
demonstrated that stress could impair the human hippocampus 
[43]. The hippocampal volume loss is demonstrated negatively 
associated with the duration of PTSD [10]. Taken together, these 
findings might imply that small hippocampal volume predisposes 
to PTSD, and PTSD in turn leads to a secondary loss of 
hippocampal volume. Although it was difficult to determine 
whether smaller hippocampal volume was a preexisting risk factor 
or a consequence of the PTSD in this study, it is plausible to 
speculate that small hippocampal volume was associated with 
chronic severe PTSD following coal mine gas explosion. 

In present study, the extent with decreased GMV in the left 
hippocampus was larger than that in the right hippocampus in 
coalminers with PTSD. The intergroup difference of the left 
hippocampal GMV was still significant after controlling the 
depression and/or anxiety, while not survive for the right 
hippocampus. Consistent with our finding, a meta-analysis found 
significant gray matter reduction of the left hippocampus in the 
PTSD patients compared with individuals exposed to trauma 
without PTSD, but no significant effect was found in the right 
hippocampus [22]. A study has shown that the left hippocampus 
was more vulnerable than the right one in the right-handed PTSD 
patients [44]. The right hippocampus appears particularly 



involved in memory for locations, whereas the left hippocampus 
is more involved in context-dependent episodic or autobiograph- 
ical memory [45]. Because re-experiencing is one of the most 
important symptoms of PTSD, the memory about the traumatic 
events seems to influence more the left than the right hippocam- 
pus. However, other studies also found greater reductions in the 
right hippocampus in PTSD patients [46,47]. The laterality of the 
hippocampal volume reduce should be further confirmed in PTSD 
patients with large sample sizes. 

In this study, no significant correlation was revealed between 
bilateral hippocampal volume and the CAPS scores, which is 
inconsistent with previous findings of the correlation between 
hippocampal volume and the symptom severities of PTSD. 
Several studies found inversely correlations between hippocampal 
volume and the measures of symptom severity [38,44,48-51], 
However, other studies had also shown no associations between 
bilateral hippocampal volume and the CAPS scores [7,9,15,18]. 
Multiple variables could attribute to these inconsistent results, 
including sample size, age and gender, trauma type and duration, 
analytic methods [19,23], and comorbid disorders [52]. The 
possible explanation for our result might be that the small 
hippocampus was a preexisting risk factor for PTSD and was only 
a trait marker but not state marker of PTSD. Other factors such as 
the small sample size and long disease course might also contribute 
to this negative result. However, further study with large sample 
size should be performed to clarify whether definitive correlation 
exists between the symptom severity and hippocampal volume in 
PTSD. 

Structural covariance between hippocampus and 
ipsilateral amygdala 

We found that the hippocampus had a deceased GMV 
covariance with the ipsilateral amygdala in coalminers with PTSD 
compared with non-traumatized coalminers. ROI-based correla- 
tion analysis revealed that this abnormal GMV covariance 
between the two regions could be attributable to the disrupted 
positive association in coalminers with PTSD 



Table 2. Hippocampal clusters showed reduced gray matter volume in PTSD patients compared with controls. 





Clusters 


Voxels 


MNI coordinates 




rvalue 


p" 






X 


y 


z 






Left hippocampus 


731 


-26 


-15 


-12 


4.390 


0.029 


Right hippocampus 


185 


23 


-24 


-12 


3.632 


0.036 



a FDR corrected at voxel level, p<0.05. MNI, Montreal Neurological Institute. 
doi:1 0.1 371 /journal.pone.01 02042.t002 
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Table 3. Intergroup comparison of GMV in bilateral hippocampal ROIs. 







Left hippocampus 


Right hippocampus 


Two-tailed f-test 


t = 4.628 


P<0.001 


f = 3.413 


P = 0.002 


ANCOVA Controlling HAMA 


F = 22.549 


P<0.001 


F= 10.617 


P = 0.002 


ANCOVA Controlling HAMD 


F = 6.954 


P = 0.012 


F=2.210 


P = 0.146 


ANCOVA Controlling HAMA and HAMD 


F- 8.081 


P = 0.007 


F=2.110 


P = 0.155 



ANCOVA, analysis of covariance; GMV, gray matter volume; HAMA, Hamilton rating scale for anxiety; HAMD, Hamilton rating scale for Depression. 
doi:1 0.1 371 /journal.pone.01 02042.t003 



Brain regions that covary with each other are often parts of a 
system that subserves particular behavioral or cognitive function 
[36]. In the healthy adults, the hippocampal volume covaries 
strongly with those of brain regions involved in the memory 
system, including the amygdala and parahippocampal, perirhinal, 
entorhinal and orbitalfrontal cortices [34] . Our study also revealed 
that the hippocampus strongly positively covaried with the 
ipsilateral amygdala in non-traumatized coalminers. However, 
the structural covariance between hippocampus and ipsilateral 
amygdala was disrupted in PTSD group. One possible explanation 
might be that there are extensive anatomical interconnections 




0.4 -I L-i — ^ 1 

L_hippocampus R_hippocampus 

Figure 2. Reduced gray matter volumes of the bilateral 
hippocampi in PTSD patients. A, the color clusters exhibit reduced 
gray matter volume of bilateral hippocampi in PTSD patients compared 
with controls (p<0.05, FDR correction). The upper numbers are the 
Montreal Neurological Institute coordinates on x axis. The colorbar 
represents the statistical value from 1 to 5. B, bar graph shows the gray 
matter volumes in bilateral hippocampal ROIs significantly reduced in 
PTSD patients compared with controls exhibiting strong effect sizes. ES, 
effect size. 

doi:1 0.1 371 /journal.pone.01 02042.g002 



between the hippocampus and amygdala [32,33]. Animal studies 
demonstrated that stress elicits hypertrophy of amygdala [53], but 
atrophy in hippocampus [54]. The structural synchronization 
between hippocampus and ipsilateral amygdala might be de- 
stroyed in PTSD. Another possible explanation for decreased 
structural covariance between hippocampus and ipsilateral amyg- 
dala was associated to the contrary functional activities of the two 
regions. Functional imaging studies on PTSD revealed increased 
activation in amygdala and decreased activation in hippocampus 
during traumatic script-driven imaginary task [55]. Growing 
dominance of amygdalar activity over the hippocampus during 
and even after chronic stress may contribute to the enhanced 
emotional symptoms seen in stress-related psychiatric disorders 
[56]. Amygdala may mediate certain stress effects on hippocampus 
[57]. The destroyed structural covariance between hippocampus 
and ipsilateral amygdala might be secondary to the long term non- 
synchronized functional activity of two regions. Overall,decreased 
GMV covariance between the hippocampus and the ipsilateral 
amygdala might be associated with abnormal processing of 
trauma-related memory such as re-experiencing, which was one 
of the most important symptoms in PTSD patients. 

Limitations 

Several limitations should be considered in this study. First, the 
relatively small sample size may impede for detecting significant 
effects because of the lack of power. However, high homogeneity 
of the PTSD patients who suffering from coal mine gas explosion 
might compensate this shortage to some extent. Second, the 
specificity of trauma event was coal mine gas explosion. Although 
patients with visible brain lesions have been excluded, we cannot 
exclude the possibility of the subtie damage caused by poisoning. 
Third, although the comorbid depression might affect the our 
results, ANCOVA result showed that intergroup GMV difference 
was still significant in the left hippocampal after controlling the 
HAMD, indicating close association between decreased left 
hippocampal GMV and PTSD. Finally, only the coalminers 
without trauma exposure were enrolled as controls. Future studies 
should also include non-PTSD controls with trauma exposure. 

Conclusion 

To our knowledge, this is the first study to investigate the 
structural covariance between the hippocampus and amygdala in 
PTSD. We found that bilateral hippocampal GMVs decreased in 
coal mine gas explosion related PTSD patients, especially the left 
hippocampus; and the hippocampal regions exhibiting reduced 
GMV had decreased structural covariance with the ipsilateral 
amygdala, which may be associated with the symptom of the 
PTSD. 
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Figure 3. Gray matter volume (GMV) covariance of hippocampal ROIs with the ipsilateral amygdala in controls. The color regions 
represent the clusters within the left (A) and right (B) amygdala in which GMV significantly co-varied with the ipsilateral hippocampal ROI (p<0.05, 
FDR corrected). The colorbar represents the statistical value. The below numbers are the Montreal Neurological Institute coordinates. 
doi:1 0.1 371 /journal.pone.01 02042.g003 




Figure 4. Inter-group differences in gray matter volume (GMV) covariance of hippocampal ROIs with the ipsilateral amygdala. The 

color regions represent the clusters within the left (A) and right (B) amygdala in which GMV covariance with the ipsilateral hippocampal ROI 
significantly decrease in PTSD patients compared with controls (p<0.05, FDR corrected). The colorbar represents the statistical value. The below 
numbers are the Montreal Neurological Institute coordinates. 
doi:1 0.1 371 /journal.pone.01 02042.g004 
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0.55 4 1 1 1 1 1 0.6 -I , , , , 

0.5 0.55 0.6 0.55 0.7 0.75 0.4 0.45 0.5 0.55 0.6 

GMV within the left hippocampal ROI GMV within the right hippocampal ROI 



Figure 5. ROI-based correlation analyses of gray matter volumes between the bilateral hippocampal and amygdala ROIs. ROI-based 
correlation analyses revealed that decreased structural covariance between the left hippocampal ROI and the left amygdala ROI (A), and between the 
right hippocampal ROI and the right amygdala ROI (B) in PTSD patients were attributed to a strong positive correlation in controls (left side: r= 0.776, 
p<0.001; right side: r = 0.568, p = 0.003) and no significant correlation in PTSD patients (left side: r = -0.205, p = 0.483; right side: r = -0.21 9, p = 0.453). 
doi:10.1371/journal.pone.0102042.g005 
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